in 



Rotation and Ha Emission Above and Below the 
Substellar Boundary 

G. Basri 

-^-, Astronomy Department, University of California, Berkeley, CA 94720 

o: 

O . Abstract. 

, I present the results of a multiyear survey of very low mass stars and 

Oh' brown dwarfs, at high spectral resolution. The spectra were gathered 

with the HIRES echelle at the Keck Observatory. Some of these objects 
are stellar and others are substellar (or ambiguous). Early indications 
that such objects can be rapidly rotating but display little Ha emis- 
sion turn out to be commonly true. This is the opposite of the relation 
^-H . between rotation and activity in solar- type stars. The Ha surface flux 

^ ' drops precipitously at the bottom of the main sequence, and seems to 

^N , be related to the luminosity or temperature of the objects. There is a 

XJ I general trend to higher rotation velocities as one looks at objects of lower 

Q>^ . luminosity. I discuss several possible explanations for these results. The 

^^ I dynamos for these objects are probably fully turbulent, driven by convec- 

^D . tion, and thus more directly related to the object's luminosity. They may 

^^ I be quenched when the rotational velocities become too fast in compar- 

,^ ■ ison to the convective velocities (supersaturation). Another possibility 

Mh| is that the atmospheres of the cooler objects are becoming sufficiently 

Q ■ neutral to decouple atmospheric motions from the field. Either of these 

^H I could explain why young brown dwarfs can be magnetically active while 

^ ■ older brown dwarfs are not. A final possibility is that instead of being 

. . I quenched, the field configuration in rapid rotators changes to a less con- 

K^ ■ ducive form for dissipative heating. This could explain why flares are 

k> I occasionally seen on generally inactive objects. 

1. Introduction 

One of the most studied properties of late-type stars is their magnetic activity. 
The quest to discover why stars of the same spectral type have such different 
levels of activity was answered in the 1980's with the "rotation-activity" con- 
nection. For solar-type stars, the faster the star spins, the greater is its total 
magnetic flux (usually measured through emission line proxies or X-ray lumi- 
nosity). The explanation for this connection is via a magnetic dynamo, which 
relies on the Coriolis force for part of its operation. The Sun is thought to have 
such a dynamo operating at the base of its convection zone (a "shell" dynamo) . 
Even for the Sun, there is clearly another dynamo mode present (Title & 
Schrijver 1998). The so-called "distributed turbulent" dynamo is thought to 
operate throughout the convection zone, and depend less on rotation. It is re- 



sponsible for most of the flux present during solar minimum. As one moves 
to later spectral types, the convection zone deepens and the size of the shell 
dynamo decreases. Thus one expects that by the time stars become fully con- 
vective, the turbulent dynamo will have taken over, and the disappearance of 
the shell dynamo is hardly noticeable. This occurs at about spectral type M3, 
and explains why activity levels do not suddenly change there (Giampapa et al 
1996). 

The field emerges on smaller scales from the turbulent dynamo than the 
shell dynamo, and are less efficient at removing angular momentum. At about 
the fully convective boundary, the spindown time of stars begins to lengthen. 
Most early M field stars are found to be slowly rotating (except when young), 
and the incidence of measurable rotation picks up at about M3.5 (Delfosse et 
al 1998). The relation between rotation and activity still appears at mid-M in 
the sense that fiare and dMe stars tend to be more rapid rotators. The fraction 
of stars which are dMe increases at later types; this is mostly because of the 
"contrast effect" whereby the same surface fiux of a given hot diagnostic stands 
out better and better against the cooling photospheric continuum. Put another 
way, a constant emission line equivalent width represents less and less surface 
flux at cooler effective temperatures. 

In the late M stars, the connection between rotation and activity is definitely 
changing. Basri and Marcy (1995) found an M9.5 star which has very rapid 
rotation, but no Ha emission. A second example was found by Basri et al 
(1996). Here we present the results of a multi-year survey of very late-type stars 
and brown dwarfs for rotation speeds and Ha emission. While most objects are 
stellar, some are confirmed brown dwarfs, and others may be substellar. The 
first part of this paper presents the results from that survey. In the second part 
we speculate a little about what the explanation for these results might be. 

2. Observational Results 

All observations were conducted with the HIRES echelle on the Keck I 10-m 
telescope. We used settings similar to those of Basri & Marcy (1995). There 
were more than 50 targets, all M5 or later, including several examples of the 
new L spectral type. Some targets were in Basri et al (1996), and a few data 
have been added from the surveys of Basri et al (1996) and Tinney and Reid 
(1998). 

2.1. Ha Emission 

Because of the extremely cool photospheres for the stars in our sample. Ha can 
only show up if there is chromospheric or coronal heating. All our objects 
later than M5 and earlier than M9 show Ha emission. The fraction of stars 
with emission is know to be high after about M5 (Hawley, Gizis & Reid 1996), 
although they find a peak fraction of 60%. Some of the emission strengths we 
measure are near or below their observational limits. At about M9.5, we see 
a marked decrease in the emission, with several stars showing no emission and 
others with only weak lines. This data is listed in Table ||, and depicted in Fig. ^. 
There is a sudden drop in the fraction of active stars at about the end of the M 
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Figure 1. Rotation, Spectral Type, and Ha emission. Each ordinate 
shows a measured rotation velocity and the abscissa gives M spectral 
subclasses. For 10 and beyond these are L spectral subclasses (plus 10). 
The size of each point is scaled logarithmically with the Ha equiva- 
lent width (small points indicate low emission). The small symbols 
represent no detected emission at all. 

sequence. The large equivalent widths are all in the late M stars (though not all 
stars are active), while the L stars are all inactive. 

A given value of emission equivalent width represents a dramatically weak- 
ening surface flux as we move into the late M and L stars. The continuum, 
which defines the normalization of equivalent width, is dropping very quickly 
with temperature (Ha now occurs in the Wien part of the Planck function). 
We have used models by Allard and Hauschildt (see Mohanty et al, this volume) 
to estimate the conversion of equivalent width to surface flux. The models in- 
dicate that the continuum flux near Ha drops by almost a factor of 2 for each 
lOOK in effective temperature. 

One of the main results of this paper is shown in Fig. g. As was pointed 
out by Basri & Marcy (1995), the fact that Ha equivalent widths remain fairly 
constant (and low) throughout this temperature range implies that the surface 
flux from heated plasma is dropping dramatically. It seems, in fact, that the 
drop in surface flux could be directly related to the temperature, or more likely 
but equivalently, the luminosity (since these dwarfs are all about the same size). 
To make it clear that the relative activity levels are falling, we show Lfja/Lboi 
in Fig. |2|. I have approximated Lboi by using our estimated temperature scale 
to spectral type conversion, and making a (small) correction for radius. With 
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Figure 2. The relation between activity and spectral type. Unlike 
with solar-type stars, they are closely tied together. Cooler stars show 
much less activity (as measured either by luminosity ratio (as here) 
or by surface flux. The symbol sizes are scaled logarithmically by 
rotation velocity. The lower boundary of points are actually upper 
limits (emission not detected). 



either surface flux or luminosity ratio, there is a dramatic decrease in the implied 
non-radiative heating (or amount of chromospheric plasma) as the objects get 
cooler and less luminous. 

One might wonder whether there is some problem with Ha in cool dusty 
atmospheres that prevents it being in emission. The point is that we are consid- 
ering chromospheric plasma, which by definition is hot enough to produce Ha . 
The absence of emission thus implies that there is no such plasma present. There 
cannot be a substantial corona in the stars showing very weak Ha emission, 
because it would create a chromosphere by photoionization (Cram 1982) that 
would easily show up. I therefore conclude that the amount of stellar activity 
drops steeply at the bottom of the main sequence, and that its levels may be 
fairly directly tied to the stellar luminosity. 

2.2. Rotational Velocities 

We determine rotational velocities by cross-correlation with rotational stan- 
dards, following the general procedures of Basri & Marcy (1995). We used 
molecular features near 8500Afor this purpose. For the M stars, GJ 406 was 
our standard. For the L stars we do not have any non-rotators in the sample, 



so we used the slowest one: 2MASS 1439+19. We checked its velocity through 
an intermediate calibration with LHS 2924. We carefully considered the effects 
of a calibrator with a finite velocity, and found that it does not substantially 
effect the determinations for v sin(i) 20 km/s or higher (slower cases may be 
overestimated, but are good to 5 km/s or better). Pressure broadening is not 
really a factor for these rapid rotators. 

Down to about M3, there are rather few field stars with v sin(i) easily 
detectable (above 3 km/s). A good summary of the situation is given by Delfosse 
et al (1998). The situation is reversed in our sample: there are rather few stars 
too slow to be measurable. And this is despite the fact that v sin(i) is always 
a lower limit on the true equatorial velocity. Our velocity determinations are 
listed in Table |l|, and shown in Fig. ||. For convenience, we refer to stars with 
V sin(i) < 10 km/s as "slow", and those with higher velocity as "fast". At M5 
in our sample, there is only 1 star which is fast. For later subclasses, the fast 
stars are as numerous as the slow ones. It is very rare to find early field M stars 
with rotations higher than 20 km/s(and they are young and active), but there 
are several later than M5 and they are not particularly active. By the end of the 
M spectral class there are no slow stars left, and that remains true for all our 
L stars (indeed, we had trouble finding one slow enough to serve as a rotational 
calibrator). 

Is the rapid rotation intrinsic, or just a result of looking at young objects? 
There is a bias against finding old stars in the L class, as there is an increasing 
fraction of brown dwarfs there. These get even cooler and fainter with age, and 
so become harder to find. Those L stars that have been shown to be substellar 
by the lithium test have all been less than 1 Gyr old. Almost none of our M 
stars show lithium, so we can assume they are mostly stellar, and older. At the 
very least, it is clear that the average rotation velocity is much higher for the 
L stars than late M stars, and the late M stars are more rapidly rotating than 
the early M stars. Spindown times obviously increase as the mass of the objects 
decreases. For the L stars, it is striking that their velocities are so high, even 
before correcting for projection effects. It seems that the typical such object 
has a rotation period measured in hours; many of them are rotating faster than 
Jupiter. 

2.3. The Rotation-Activity Connection 

For solar- type stars, a clear increase in magnetic activity is seen with an increase 
in rotation velocity. A nice summary figure for the rotation-activity connection 
can be found in Randich 1998 (see also Saar, this volume). The variables in 
which this relation is best expressed are either surface flux vs. rotation period, 
or luminosity ratio (with \-iboi) vs. Rossby number. The Rossby number is itself 
the ratio of convective overturn time to rotation period. Fig. 11] shows that the 
activity levels do not increase with rotation velocity. There seems to be little 
connection between the two at low velocities, and none of the rapid rotators 
have high emission levels. This is without regard to the spectral type. 

The dependence on spectral type can be seen in Figs. || and ^. They show 
that the L stars are both rapidly rotating and inactive, while late M stars span 
a range of rotations but are more active. The L star sample is substantially 
smaller, so it is possible that more slowly rotating or more active examples will 



Table 1. Observed Rotation and Ha Emission. 



Star 


Spectral Type 


V sin(i) Ha Eq.W. 






(km/s) 


(A) 


GJ 1093 


M5.0 


3.0 


1.0 


GJ 1230B 


M5.0 


<3.0 


<0.2 


GJ 1057 


M5.0 


<3.0 


<0.2 


GJ 905 


M5.0 


<3.0 


-0.2 


GJ 810B 


M5.0 


<3.0 


-0.6 


GJ 3076 


M5.0 


15. 


7.0 


GJ 3454 


M5.0 


8.0 


7.0 


GJ 1156 


M5.0 


9.0 


6.9 


GJ 1154A 


M5.0 


5.0 


5.9 


GJ 273 


M5.0 


<3.0 


-0.2 


GJ 65AB 


M5.5 


30. 


4.4 


LP 759-25 


M5.5 


12.5 


2.5 


GJ 1245A 


M5.5 


17.5 


4.0 


GJ 1245B 


M5.5 


7.0 


4.4 


GJ 1002 


M5.5 


<3.0 


<0.2 


GJ 1245A 


M5.5 


25. 


3.5 


YZ CMi 


M5.5 


7.0 


6.5 


GJ 3828B 


M6.0 


6.0 


1.1 


GJ 412B 


M6.0 


8.0 


9.4 


GJ406 


M6.0 


<3.0 


5.3 


GJ 1111 


M6.0 


13. 


4.3 


GJ 3622 


M6.0 


<3.0 


2.9 


GJ 2005 


M6.0 


9.0 


3.0 


CTI 2332+27 


M6.0 


25. 


6.3 


LP 713-47 


M6.5 


13. 


4.2 


CTI 1747+28 


M6.5 


45. 


5.0 


CTI 0156+28 


M6.5 


5.0 


2.3 


CTI 1539+28 


M6.5 


8.0 


0.5 


GJ 4281 


M6.5 


6.0 


4.3 


VB 8 


M7.0 


5.5 


5.5 


CTI 1156+28 


M7.0 


11.5 


2.8 


GJ 3877 


M7.0 


<3.0 


4.0 


LHS 2645 


M7.5 


6.5 


4.9 


LHS 2632 


M7.5 


<3.0 


0.4 


2MASS 1254+25 


M7.5 


15. 


7.2 


2MASS 1256+28 


M7.5 


7.0 


3.0 


LP 412-31 


M8.0 


9.0 


18.4 


RG 0050-2722 


M8.0 


<5.0 


2.9 


LHS 2243 


M8.0 


<5.0 


15.8 


VB 10 


M8.0 


8.0 


4.3 


LHS 2397a 


M8.0 


22.5 


20.5 


2MASS 1242+29 


M8.0 


5.0 


8.1 


TVLM 513-46546 


M8.5 


60. 


1.8 


CTI 1156+38 


M8.5 


11. 


3.7 


DENIS 2146-21 


M9.0 


10. 


<0.2 


LHS 2924 


M9.0 


10. 


1.5 


CTI 0126+57 


M9.0 


10.5 


4.2 


TVLM 868-110639 


M9.0 


37.5 


9.6 


DENIS 1208+01 


M9.0 


10. 


- 


LHS 2065 


M9.0 


9.0 


8.4 


BRI 1222-1222 


M9.0 


2.0 


9.7 


ESQ 207-61 


M9.5 


10. 


4.7 


DENIS 0021-42 


M9.5 


12.5 


0.5 


LP 944-20 M9.5 


M9.5 


28. 


1.0 


BRI 0021-0214 


M9.5 


37.5 


<0.2 


DENIS 0909-06 


LO. 


25. 


<0.2 


G196-3B 


LI. 


10. 


<0.2 


2MASS 1439+19 


LI. 


10. 


<0.3 


Kelu-1 


L2. 


60. 


1.5 


GD 165B 


L3. 


37.5 


<0.2 


DENIS 1058-15 


L3. 


32.5 


<0.3 


2MASS 1146+22 


L4. 


22. 


<0.2 


DENIS 1228-15 


L4. 


11. 


0.5 


LHS 102B 


L4. 


32.5 


1.6 


DENIS 0205-11 


L5. 


22. 


<0.2 


2MASS 1632+19 


L5.5 


30. 


<0.2 


DENIS 0255-47 


L6.D 


40. 


<0.2 



be uncovered. But it is already clear that they are different. Based on the 
lithium test, at least a third of them are brown dwarfs. This fraction increases 
toward later L types, and becomes 100% by about L4 or L5 (where we think the 
minimum main sequence temperature lies). 

The general conclusion that can be drawn from this is that the rotation- 
activity connection is weakened once one crosses the fully convective boundary, 
and that it completely disappears (and perhaps even reverses) cooler than M9. 
Since spindown is due to angular momentum loss via magnetic winds, it is not 
surprising that low activity levels go along with rapid rotation. To the extent 
that some of our objects are on the main sequence, they are probably old and 
show that there is little spindown over long times. In order to test older brown 
dwarfs, we will have to measure rotation velocities for the methane dwarfs. 



3. Possible Explanations for the Results 

There are several possible explanations for these results. The truth may in- 
volve more than one of them (or another not discussed here). One is that the 
ionization levels in the photosphere may have become so low that there is in- 
sufficient conductivity to allow coupling of the magnetic field to the gas. Then 
gas motions may not twist up or excite waves in the fields, and there is no dis- 
sipation to heat the upper atmosphere. This has to be true even in the face of 
ambipolar diffusion, which couples small numbers of ions to the neutrals fairly 
effectively (as in T Tauri disks). The alkali metals, which are the last source of 
electrons, are becoming quite neutral in the L stars. A possible counter-example 
to this hypothesis is provided by the detection of (non-flaring) Ha emission in 
a methane dwarf, which should be very neutral (Liebert et al, this volume). 
Another potential wrinkle is that the plasma is hotter beneath the photosphere 
and will become ionized at some depth. There will be convective motions there 
too, so it is not obvious how the surface fields will behave. 

The traditional rotation-activity connection may arise because activity in- 
creases with decreasing Rossby number. In solar-type stars, activity levels in- 
creases steadily from a Rossby number of 5 down to 0.1. Coronal activity "sat- 
urates" for Rossby numbers between about 0.1-0.01; the activity flattens out 
at its maximal values seen. It has been shown that most activity diagnostics 
scale simply like L^/ Lfeoj (Schrijver 1987). All low mass objects should have 
turbulent dynamos, which are driven by convective motions. Rotation can en- 
hance production of fields, and the amplitude of convective velocities also does. 
But convective overturn times scale with effective temperature in these objects. 
At the bottom of the main sequence they can increase to months, while typical 
spin periods are dropping to hours. A simplistic analysis using mixing length 
theory suggests that the Rossby numbers for stars in our sample lie mostly in 
the range 0.01-0.001. This is because the convective overturn times scale down 
with effective temperature, and the rotation velocities get higher for the cooler 
objects. 

As first noted by Vilhu & Rucinski (1983) and nicely summarized by Randich 
1998, there is some evidence for "supersaturation" - an actual turndown of ac- 
tivity levels - for Rossby numbers less than 0.01. I speculate that the dynamo 
may be unable to operate efficiently at such small Rossby numbers, perhaps be- 



cause rotation organizes the flows too much. A possible counterexample to this 
hypothesis is provided by the very rapid rotator Kelu-1, which exhibits a per- 
sistent (though rather weak) Ha emission line. This hypothesis can be tested 
when a sufficient sample of L stars with Ha emission has been collected. They 
should show systematically lower rotational velocities if rotation is important in 
this way. 

The only population of brown dwarfs that is known to be fairly active is 
the very young ones. Those found in star forming regions (eg. Wilking et al 
1999) or as old as the Pleiades (eg. Zapatero-Osorio et al 1997) tend to show 
reasonably strong Ha , or even X-ray emission (Neuhauser et al 1998). These 
are also of course more luminous, and perhaps their convective velocities are still 
high enough to prevent supersaturation. The Pleiades objects have been shown 
to be generally rapid rotators (Oppenheimer et al 1998). 

A related possibility is that the field is not actually quenched by rapid 
rotation, but instead takes on a relatively stable, low multipole character like 
that of Jupiter. In that case, the field might be sufficiently quiet (especially 
in conjunction with low conductivity) that it does not suffer the dissipative 
configurations that power stellar activity. Thus, the objects might still have 
strong fields, but little stellar activity. This could be tested in principle using 
Zeeman diagnostics. Valenti (this volume) suggests using FeH for objects in 
this temperature range, and have shown that it can work in late M stars. A 
suggestion that there are indeed fields comes from the fact that flares have been 
seen in objects that seem otherwise quiescent, such as VBIO (Linsky et al 1995) 
and 2MASSW J1145572+231730 (Leibert et al 1999). 

It is clear that a new regime in the study of stellar activity has been reached 
at the bottom of the main sequence. This will be the motivation for much 
additional work in the coming years as our sample of very cool objects is greatly 
increased by the 2MASS, DENIS, and Sloan sky surveys. We will have to follow 
up with low and high dispersion spectroscopy and synoptic photometry to gather 
the data needed to fully sort this problem out. But in the process, we are likely 
to gain valuable new information on both stellar dynamos and the process of 
non-radiative heating that produces stellar activity. 
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Discussion 

C. Bailer- J ones: I am a bit concerned about your v sin(i) measurements. How 
do you deduce v sin(i) in the absence of a well-known slow enough and late 
enough standard? 

G. Basri: (This is commented on in the manuscript) We worried about this 
quite a bit. Originally we developed a means of correcting the answer using a 



moderately fast rotator. But then we found one at 10 km/s which is plenty slow 
enough for most of the L stars. 

R. Jefferies: Are the phenomena you have presented related to the question of 
the period gap in CVs? 

G. Basri: Yes, they probably are. The same sort of conductivity arguments 
have been invoked recently by the Meyers for the period gap systems, and the 
low mass secondaries are presumably similar to the objects I discussed. Indeed, 
the rotation periods are sometimes near 2 hours as well. 
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